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ZAjhe  bum  stability  of  a  thermonuclear  reacting  plasma  is  examined  under  the  assumption  that 
a  bum  equilibrium  exists  due  to  the  rapid  increase  of  loss  rate  with  plasma  beta  once  a  critical 
beta  value  is  exceeded.  It  is  found  that  perturbations  about  equilibrium  generally  result  in  a 
rapidly  damped  exponential  decay  and  a  second  slow  root  which  can  be  either  growing  (unstable) 
or  damped.  In  the  case  where  the  slow  root  is  growing,  the  possibility  that  it  can  be  stabilized 
by  feedback  control  of  the  rate  at  which  neutral  gas  Is  fed  into  the  system  is  considered. 
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STABILITY  OF  BETA  LIMITED  THERMONUCLEAR  BURN 


Plasma  heaClng  by  fusion  products  Increases  with  temperature. 
Thus,  once  the  Ignition  point  Is  exceeded.  In  order  to  establish  a 
thermonuclear  burn  equilibrium,  it  is  necessary  that  plasma  losses 
increase  sufficiently  rapidly  with  temperature.  In  this  note  we 
examine  the  implications  of  the  hypothesis  that  such  an  equilibrium 
exists  due  to  a  rapid  increase  in  losses  with  plasma  beta  once  a 
critical  value  of  beta  is  exceeded.  For  example,  in  tokamaks  such  a 
critical  beta  might  result  due  to  the  onset  of  magnetohydrodynamlc 
instabilities.  Since  generally  accepted  theoretical  models  for  anoma¬ 
lous  plasma  loss  in  such  devices  as  tokamaks  are  not  available,  we 
shall  try  to  make  as  few  assumptions  as  possible  regarding  the  depend¬ 
ences  of  loss  rates  on  plasma  parameters  except  for  the  existence  of 
a  critical  beta. 

We  utilize  the  following  zero  dimensional  model  equations  for 
a  deuterium-tritium  mixture. 


N0  -  -  Na  +  8NdNt  , 


8D  "  -  ^N  ND  +  SD  -  %NT  * 

nt  "  -  yJ  nt  +  sT  -  0Ndnt  , 


Note:  Manuscript  submitted  November  21,  1979. 
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B  -  -  Y*NeT  -  Y®NdT  -  Yj»tT  -  yX*  +  6E„Vl  •  <U) 

where  is  the  particle  number  of  species  o  (o  ■  a,  D,  T,  e  for  alpha 

particles,  deuterons,  tritons,  and  electrons),  Y»  is  the  loss  rate  of 

N 

particles  for  species  o,  yf  is  the  energy  loss  rate  for  species  o,  8(T) 

Et 

is  the  deuterium-tritium  thermonuclear  rate  coefficient,  E  is  the 

a 

alpha  particle  creation  energy,  E  ■  (Ng  +  +  N  JT  is  the  thermal 

particle  energy  content  of  the  plasma,  T  is  the  temperature  which  has 
been  taken  to  be  equal  for  all  particles  (this  is  approximately  true  if 
the  equilibration  times  are  sufficiently  short  compared  to  the  confine¬ 
ment  time),  and  the  dot  denotes  time  derivitive.  Ohmic  heating  has 
been  neglected  in  Eq.  (Id)  since  its  contribution  is  small  at  the  burn 
point.  The  quantity  S0  denotes  the  rate  at  which  particles  are 
supplied  from  an  external  neutral  gas  input  (the  effect  of  recycling 
of  particles  transported  to  the  edge  may  be  absorbed  in  the  coefficient 


D  T 

Assuming  yn  »  yn  and  SD  *  ST»  Eqs.  (lb)  and  (lc)  imply  that 
(Np  -  N^,)  'v  exp  (“JN^  <&)•  Thus  under  these  circumstances  it  is  rea¬ 
sonable  to  take  Np  »  Nj  and  eliminate  Eq.  (lc).  Assuming  Na  <<  N^.Np, 
neglecting  the  term  on  the  right  hand  side  of  (lb),  and  utilising 

the  neutrality  condition,  Ng  ■  +  2Na,  we  obtain  from  Eqs.  (1) 


ND  "  -  +  SD  * 


T  +  &eahD  , 


In  addition,  we 


where  ea  h  Ea/4  and  Ye  =  \  yJ  +  J  Y J  +  J  yJ  -  yJ. 
shall  allow  for  feedback  control  of  the  external  gaa  puffing  rate,  Sp 
We  assume  that,  based  upon  sensing  of  the  particle  number,  Np,  the 
gas  supply  rate  is  adjusted, 

SD  -  -  v  [Sp  -  PCNp)]  ,  *  (2c) 

where  v-1  is  a  time  constant.  This  time  constant  reflects  the  fact 
that  the  control  system  will  have  a  finite  response  time,  as  well  as 
the  fact  that  a  particle  supplied  at  the  edge  will  take  a  finite  time 
to  diffuse  into  the  reacting  core  plasma  region.  The  latter  time  is 
probably  of  the  order  of  ^y^  and  sets  a  upper  limit  on  v.  For 
simplicity  we  shall,  in  what  follows,  neglect  any  T  or  Up  dependence 
of  v. 

From  Eqs.  (2a)  and  (2b)  the  equilibrium  burn  condition  is 

s„  -  rj  V  (3.) 

?  »D  *  Ge  +  VV*  •  C3b) 

where  the  overbar  denotes  equilibrium.  Linearizing  Eqs.  (2)  for 
perturbations  about  this  equilibrium  of  the  form  exp  (pY^t)  we  obtain 

(p  +  1  +  tp  +  (y  +  1)  (1  -  Oj)  +  Yo4)  (P  +  v)  + 


a,a,v  -  a,(p  +  v)  (Ya,  -  Y  -  2)  +  a,v[p  +  (Y  +  1)  (1  -  «=)+  Y<*,]  -  0, 


where 

Y  "  Ye/YN»  V  “  ^/yN*  °1  "  (5D/YH)  3yN/3ND»  °2  "  <5/yS>  3yN/yT’ 
«3  -  (%/ve)  3Ye/3ND,  «4  -  (T/Ye)  3Ye/8T,  atj  -  (T/S)dB(T)/dT,  end 

a,  -  (N_/S_)  dF(Nn)/dNn.  In  the  case  where  there  is  no  feedback  and  S 
is  held  constant  independent  of  plasma  parameters,  Eq.(4)  becomes  (v*0) 

(p  +  1  +  a^)  [p  +  (y  +  1)  (1  -  «5)  +  Y<*4J  -  «2<Y a3  -  Y  -  2)  -  0. 

(5) 

We  now  consider  a  special  case:  there  is  a  rapid  Increase  of 
energy  loss  (y£)  past  the  critical  value  of  B»  but  the  particle  loss 
rate  (yj!)  is  relatively  uneffected  (this  situation  might  result  from 
the  onset  of  magnetic  field  braiding  in  a  tokamak). 

In  this  case 

a3  %  a4  >3>  °1  ^  a2  ^  1  (6) 

(Note  that  since  the  a  terms  are  derivative  of  logarithms  with 
respect  to  logarithms,  much  greater  than  inequalities,  such  as  (6), 
are  in  actual  situations  likely  to  be  marginal  at  best.  Thus  the 
case  which  we  consider  should  be  regarded  as  illustrative  only.) 

For  this  situation,  Eqs.  (5)  and  (6)  show  that  without 
feedback  there  are  two  roots  p  -  p^,  p2. 


P i  *  -  y«4  , 

(7a) 

P2  (®2flf3^ar4^  "  ^  + 

(7b) 
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Th«  root  represents  a  rapidly  daapad  solution  (-p^  »  |  p2  |  ) . 
Ths  root  p2  can  ba  aithar  daapad  or  growing  dapandlng  on  tha  signs 
of  a2  and  and  tha  ralatlva  nagnltudas  of  tha  a  coafficlants.  We 
now  consider  tha  possibility  of  stabilising  tha  root  p2  by  feedback 
(for  the  case  p2  >  0).  Applying  (6)  and  (4)  we  find  that  stability 
results  if  both 

v  >  p2  and  dT/dH^  >  P2  ,  (8) 

and  the  burn  is  unstable  otherwise.  Thus  we  see  that  stabilisation 
occurs  if  the  system  response  tlae  is  less  than  the  instability 
growth  tlae,  and  If,  in  addition,  the  feedback  Is  strong  enough 
(dF/dNp  >  p2y^).  Since  a  lower  llait  on  y  results  from  the  finite 
time  for  a  particle  supplied  at  the  edge  to  reach  the  reacting  core, 
it  is  not  clear  that  (8)  can  be  satisfied. 
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